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In telephone traffic studies, the observed proportion of unsuccessful at- 
tempts over a given time interval is one of the measures commonly used to 
evaluate the grade of service provided by trunk groups. This paper deals 
with the derivation of an approximate formula for the variance of this 
estimate when (i) call arrivals constitute a Poisson process, (ii) service 
times are independent of each other and identically distributed according 
to a negative exponential law, and (Hi) calls placed when all trunks are 
busy are either cancelled or sent via some alternate route (loss system). 
Comparison of simulation data with numerical values computed by means 
of this formula indicates that the latter is accurate enough for practical 
purposes. 

The observed proportion of lime during which all trunks are occupied is 
also an estimate of the grade of service (defined as the probability that a 
call will be lost or overflow). It is shown here, that for relatively small loads, 
this estimate has a smaller variance than the observed proportion of lost or 
rerouted calls. However, as the load is increased, the inequality between the 
variances of these two estimates is reversed, the cross-over occurring in the 
vicinity of the point where the load (in erlangs) is equal to the number of 

trunks. 

For a given observation period, the proportion of time when all trunks are 
busy can be either measured exactly or estimated by "switch-counting." 
In the latter case, the group is scanned at regidar intervals and one ob- 
serves, for each scan, whether all trunks are occupied or not. The average 
number of scans which indicate that all trunks are busy is an estimate of 
this proportion and is, a fortiori, an estimate of the probability of loss. 
The effect of the scanning rate on the accuracy of this estimate is investi- 
gated. 

I. INTRODUCTION 

In this paper, we shall consider the simplest type of loss systems, 
namely full availability groups with Poisson inputs, negative exponential 
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service times, and cancellation or rerouting of calls finding all trunks 
occupied. Under these assumptions, we shall obtain an approximate 
expression for the variance of the measured call congestion, the latter 
being defined here as the proportion of calls which either are lost or 
overflow to some alternate group during a given time interval. In the 
derivation of this expression, use is made of the classical formula for the 
propagation of errors, whose computation requires the evaluation of 
the first- and second-order moments of the joint distribution of the 
number of offered and the number of overflow calls. Since the marginal 
means and variances of this distribution are known (cf. Ref. 1), the 
emphasis is placed here on the derivation of the covariance. Computed 
values of the variance of the measured call congestion are shown to be 
in good agreement with simulation results (cf. Figs. 1-5).* Charts 
giving the variance of this ratio for group sizes up to 50 and offered 
loads (in erlangs) per trunk of 0.1 to 10, are reproduced in Figs. 6-8.* 

For a given observation period, the measured call congestion and the 
observed proportion of time when all servers are busy — here called 
measured time congestion — provide us with two estimates of the 
probability that a call will either be lost or overflow to some alternate 
route. Neither of these two estimates has a uniformly smaller variance 
than the other. Actually, the following holds: for relatively small loads, 
the measured time congestion has a smaller variance than the measured 
call congestion. However, as the load is increased, the direction of the 
inequality is reversed, the cross-over occurring in the vicinity of the 
point where the load (in erlangs) is equal to the number of trunks. Thus, 
the measured time congestion is not always a more efficient estimator 
of the probability of loss than the measured call congestion (cf. Figs. 
9 and 10). 

(In what follows, the terms measured time congestion and measured 
call congestion will always be abbreviated to time and call congestion, 
respectively. These terms will refer throughout to measurements per- 
formed over a given time interval.) 

For a given observation period, time congestion can be either meas- 
ured exactly or estimated by switchcounting. In the latter case, the 
group is scanned at regular intervals and one observes, for each scan, 
whether all trunks are busy or not. The proportion of scans which indi- 
cate that all trunks are busy is an unbiased estimate of time congestion 
and is, a fortiori, an estimate of the probability of loss. Clearly, the vari- 
ance of this estimate increases as the scanning rate decreases. The loss 
of accuracy due to scanning is depicted in Fig. 11. 

* See illustrations placed later in this article. 
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Under the present assumptions, loss probabilities can also be esti- 
mated from carried loads measured either exactly or by scanning. For 
offered loads (in erlangs) falling short of the number of trunks, simula- 
tion has shown that such estimates have smaller variances than the 
estimates mentioned earlier. This fact is illustrated in Fig. 12. The 
effect of scanning on the accuracy of loss estimates based on load 
measurements is sketched in the same figure. 

Finally, we note that estimates of loss probabilities based either on 
observed call congestion or on carried load measurements are biased, 
respectively, downwards and upwards. These biases are, however, 
quite small and likely to be negligible in most situations of practical 
interest. 

II. THE COVARIANCE FUNCTION 

Consider a group of c trunks which operate in parallel and are fully 
available to all requests. If a call is placed when a trunk is free, service 
starts immediately; otherwise the request is either cancelled or routed 
via some alternate group (loss system). Regarding the input and the 
service durations, the following assumptions will be made : 

(i) The time intervals between successive service demands (whether 
successful or not) are independent of each other and have a common 
negative exponential distribution with mean equal to \/a (Poisson 
input). 

(it) The service times are independent of each other and have a 
common negative exponential distribution whose mean will be taken 
throughout as the unit of time (a is therefore the offered load in erlangs). 

The following notation will be used: 

N(t) = number of busy trunks at time t, 

R(t) = total number of requests offered during (0,0 , 

S(t) = total number of unsuccessful requests during 

(0,0, 

P(t,n,r,s) = Pr [N(t) = n, R(t) = r, S(t) = s]. 

From the definition of P (-,-,•,• ), it follows that: 

P(t,n,r,s) = for n > c, (t ^ 0) 

P(t,n,r,s) = Ofors > r, (t ^ 0) 

P(Q,n,r,s) = Oforr ^ 1. 

It will be convenient to extend the definition of P (-,-,•,• ) and to 
adopt the convention : 
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pansion of F\ in powers of x and the summation on the right of (6) 
is the coefficient of x c in that same expansion. Therefore 

G( ) . _ (« + fajM - gW (7) 

(c + 1)(7 C+1 (w,j/) - ayz-<r c (w,y) 
where the t's and o-'s are defined by 

e avx (l-xy lv - 1) -«= £>.(«*)*" (8) 

o 

Ke ix ~ 1)av (l - x ) a(v - 1} - w f ' u a{1 - v)+w - l e a{y - l)u du = j^r n (w,y)x n . (9) 

■'0 

We note that if c n (-,a) and L n (a) (•) stand respectively for the Pois- 
son-Charlier and the Laguerre polynomials of degree n and parameter 
a; i.e., if (Ref. 2, pp. 34-35 and 101, and Ref. 3, p. 26) 



n (t,a) = «»">!)* g (-l) n ^y-'t(t - 1) •• 

md 

L, M (l) =±( n + a )tf. 
t=h \n + v ) v \ 



(t ~ v + 1) 



and 



then: 

e ax (l -x) 1 = itcn(t,a)[(ax) n /?ll] 



= 2 (-i)"C(«K 



*n(«w) = [(aZ/)Vw!]c„[a(?/ - 1) - w,a?/] 

= (-l)"L„ [a(2/ - 1) -'°- n] (a2/). 

[The relation between the o-'s and Kosten's ^-functions (cf. Ref. 1) is 
readily found. Indeed, by definition 

00 

a(x — 1) /i \z V - * z n 

e (1 — Z) = 2-HPnX 



so that 

< r B («;, 2 /)= e V <l(J '- 1) -' .] 
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(10) 



For later purposes we note that: 

<r m+1 (w,y) = am+i(w + l,y) — <r m (w + l,?/), 

(m = 0, 1, •••) 
[w - a(y - l)]«r m (w + 1,0 ) = (m + l)a m+ i(w,y) 

- aya m (w,y), (11) 

(m = 0, 1, •••) 

c 

S <Tn(w,y) = <T C (W + 1,?/). (12) 



These identities are immediate consequences of recurrence relations 
known to hold for the corresponding Laguerre polynomials (cf. Ref. 2, 
p. 98). 

Substituting (7), (8) and (9) into (5) yields: 

i? f \ (c + 1) Tc+i (w,y ) - ayz-T c (w,y) v , * „ 

(c + l)a c+ i(w,y) - ayz-a c (w,y) ^ (13) 

+ £ T n (w,y)x". 

Wc can now obtain the generating function, F (■,-,■), of the Laplace 
transforms of the joint probabilities Pr [R(t) = r, S(t) = s], r,s = 0, 
1, • • • , by deleting from (12) all terms of degree higher than c in x 
and then setting x equal to 1. If we perforin these operations and then 
make use of (12), wc find that: 

n , n (c + l)T c+i (w,y ) - ayz-T c (w,y) , . 

(c + 1)<T C+I (w,y) - ayz-a c (iv,y) 

(14) 

+ X) T n (w,y) • 



The moments of the joint distribution of R and S can now be obtained 
by evaluating the derivatives of F(w,',-) for y =■ z — 1. 

Differentiating (7) with respect to z and making use of (11), we find 

ay-T c (w,y) 






=i w — a(y — 1) 

_ ay ■ (T c (w,y) [(c + 1 ) t c+1 (w,y) — ay-r c (w,y) ] 
[w - a(y - 1)1 2 <t c (w + \,y) 

In particular, for y = 1 , we have the well-known result 

dF _ aE i, c (a) 

dz „=*=i w 2 



(15) 
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so that 

ES(t) = atE lt „(a) (16) 

where Ei iC (a) is Erlang's loss formula. 

Taking the derivative of (15) with respect to y and then setting y 
equal to 1, yields: 



d 2 F 



dydz 



= - T C (W) 1 + - + -. r-rr 

W Wa c {W 4- 1)J 



=*=i W 



wdy 



a<r c (w) [, . .\ fl , N 

1) 5T. r «+lW) 



(17) 



w 2 <t c (w + i) L °y 

- a — r c (w,y) 
dy J 



where 

a m (w) = a,„(w,l) and r m (w) = T m (w,l), (m = 0, 1, • • • ). 

To determine the derivatives of T c (w,y) and r c +\(w,y) with respect 
to ?/, consider the generating function 

H(w,x,y) ==Ke"- 1)a " (1 - .r) a(w " ,)_ '" / " u a *~ v) * B ~ 1 e a(v ~ 1)u du. 

•'n 

Differentiating this expression with respect to y and then setting y 
equal to 1, we find that 

— H(w,x,y) \ v =i= 5Z j- r m (w,y)x n \ y=x 
ay o d# 

a(i-l) 

= ^--^H— t(1+««) 
w 2 (l + w>) 

and, therefore: 

I r m (w,y) \ V ^ = K {m _ 1 ~ w(1 + w) (l + ^) » (lg) 

(ro = 1, •••)• 

Since 

Xe"V 



T m 



(«) = 



w-m! 
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we obtain, upon taking (18) into account: 

£|£[R«) -flCOll = 2K ^ + K e -^ • 2H«g> (19) 

c!w 3 c!u> 2 a w+ i(c) 

where the notation £{/| is used to designate the Laplace transform of /. 
Since ER(t) - at and ES(t) = at E iiC (a), we also have: 

*0~M,B m -Ki£.%±* (20) 

where Gov \R(t),S(t)] stands for the co variance between R(t) and S(t). 
For y = 1, m = c and w replaced by w + 1, (11) reduces to 

(w + l)<r e (to + 2) = (c + l)a e+ i(w + 1) - flkr c (w + 1), 

and (20) can be rewritten as follows: 

£[Cov [R(t), S(t)]} 

-a c+l (21) 

= K T-T7 f iw -l. ^ [(C + D*efi(w+ 1) - a*e(w + l)]. 

C\W 2 {W + 1)(T C (W) + 1) 

Let Wi , i = 1, • • • , c be the c roots of a c (w + 1). It is well known 
that these roots are simple, smaller than —1 and at least one unit 
apart. Then expanding (21) in partial fractions and making use of the 
relation (c + l)o- c+ i(0) - aa c (0) = 0, which is (11) for y - 1, W = 0, 
we find: 

£{Cov [R(t),S(t)]\ 



= K 



» o+i r 



c! L w 2 <r c (l) 

C+l d <T C+1 (W + 1) 



W dw a c (w + 1) 



_ (c + lK+xd) - q<r e (l) 

u>=0 l«0- c (l) 



+ ( c + i) . y go+xCtti+jj ] 

and the covariance between R and S is, therefore, given by 
Gov [R (0,5(0] 

= g e ~° aC+1 f (c ± ^^id) - Prod) . f 



H j 



c! L "■•(!) 

( C +l)<r c+1 (l)-acr c (l) ( 22 ) 



+ fr + D *»-h(«H i' 



dw> <t c (w> + 1) 



To(l) 



+ ( c + l)iV ^ l(tti+1)g '" 1. 
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To determine explicitly the derivative appearing in (22), let us con- 
sider Cov [R(t),S(t)] for small values of t. Writing P(t,r,s) for the 
(equilibrium) probability that, during a time interval of length t, r 
requests arrived and that, among these r requests, s of them found all 
the trunks busy, we have (t small): 

P(t,0,0) - 1 - at + o(t) 

P(t,l,0) = at [I - E lte (a)] + o(t) 

P(M,D = atEiAa) + o(t) 

P(tfl,8) = 0, (s^l) 

P(t,r,0) = o(t), (r > 1) 

and 

00 00 00 

^ T. sP(t,r,s) < 2 rP(t,r,s) < £ rsP(t,r,s) 

r,f-2 r,«=2 r,a=2 

- V 2 —at \0t) „ f.\ 

<2,re — r- = o{t). 

Hence Cov [R(t),S(t)] = at 7?,, c (a) + o(t). Letting t tend to in (22), 
we find that 

/ . 1N d ffc+i(w + 1) 



dw a c {w + 1) 

(c+ l)o- e+ i(l) - og e (l) 
••.(1) 

_ (e + DlT ac+l(wi + 1} . 

v ^ ; ' ^ w*(l + Wi) II (w»- - «»i) 

Substituting this expression in (22) yields 



Cov [R(t),S(t)] = K 



ft 



e" 1 a c+1 T (c + l)cr c+1 (l) - aa e (l) 



_ r c+ i)iy **i(«< + D wn 

^ i; • Yw/d + wo n (w* - «,-) < 22 > 
^ KC ^ XJ - v^(i + ^o n («»* - toy) J 



ACCURACY OF LOSS ESTIMATES 1149 

We shall now determine the constant term appearing in (22'). To this 
end, we note that 



1 c 
,{w + 1) = -7 TT (w — Wi) 



ci r=i 

and, therefore 

<r c+l (w+ 1) = 1 (r c+ i(l) _ 1 <r c+1 (0) 

w(l + w)a c (w + 1) u> CT C (1) 1 + w <r c (0) 

_j_ e |V 1 (Tc+l (tgj + 1) 

■ i=i(w — Wi) Wi(i + Wi) n ( w < — w i) 

Hence, for w = 0, we have: 

fc + D!± ... '-"^^ 5 

<=i u>,-(l + «;,-) II (u>i — ^i) 

= (e + 1) lim I" 1 - ^ - * SB® (23) 

«-0 \_W <T C (1) 1 + W O- e (0) 



CT c+ l(w> + 1) 



»] 



1^(1 + w)0- c (w; + 



= — a. 



Furthermore, forw = 1 + w, ; , ?/ = 1 and m = c, (11) yields: 

(Wi + IKK + 2) = (c + l)*e+l(Wi + 1). (24) 

We also note that the c roots of <r c (w -\- 2) are w { — I, i = 1, ■ ■ • , 
c, so that 

c e (w 4-2) = L f[ (w ~ Wi + 1). (25) 

Hence, combining (24) and (25), we have: 

(c + 1) Wo, + 1) _ ft _ 

1 + w.- y=i 

Using (23), (26) and the relations 

(c + 1K +1 (1) = <r e (l){c + 1 + a^i.e(o)} 



(26) 
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(22' ) can be simplified as follows : 

Cov [R(t),S(t)] = aBiM I {c -HI — o[l - E Uc {a)\)t + a 

(27) 

-§£n( 1+ ^)]- 



Since, as pointed out above: 



max Wi < — 1 

1 < i < c 



we have : 



and 



mm I Wi — Wj I > 1 

lgt'.igc 



C Uljt / -I \ 

<=i w^ f£ \ Wi - wj/ 



dt H w? jA V wi - Wi! 



c w t t / 

Hence we have the following inequalities [use is also made here of (23)]: 
aE lie (a){c + 1 - a[l - #i..(a)]}* 

< Cov [R (t),S (t)] 

< atfi.e(a) [{c + 1 - o[l - ^i..(a)]}< + a] 

and, for large values of £(>0): 
Cov [R(t),S(t)] 

= aE Uc (a) [[c + 1 - o[l - E liC (a)}}t + a] + o(e - '). 

III. VARIANCE OF CALL CONGESTION 

In the preceding section, exact and asymptotic formulas were obtained 
for the covariance between the number of offered calls, R(t), and the 
number of overflow calls, S(t), during a time interval of length t. These 
expressions can now be combined with known formulas for the means 
and variances of R(t) and S(t) to obtain an approximate expression for 
the variance of call congestion. Indeed, according to the classical formula 
for the propagation of errors, we have: 



ACCURACY OF LOSS ESTIMATES 1151 

Var [S (t)/R(t)] 



~{l/ER(t)f Var [S(t)} 4- \ES(t)/[ER(t)] 2 } 2 Var [R(t)] 

- 2\ES(t)/[ER(t)} 3 } Cov [R(t),S(t)) 
= (\/atf Var [S(t)] + E lt 2 {a)/(at) 

- 2[Ei Aa)/ (at) 2 ] Cov [R(t),S(t)] 
< (l/at) 2 Yar[S(t)] 

where, assuming t large: 

Var [S(t)] ~ at E llC {a) |~1 + 2a ~ * fl( " } n 1 

L dw o- (w 4- 1) Jw-0 

and 



(29) 



(30) 



Cov [R(t),S(t)] ~ at Ei, e (a)[c + 1 - a + aEi..(a)]. (31) 

We note that Var [S (i)//2(<)] is, asymptotically, of the form fc/J, 
where k depends only on a and c. 

The exact and asymptotic expressions for Var [S(t)] were first de- 
rived by Kosten, Manning and Garwood. 1 These formulas can be ob- 
tained in a straightforward manner from the generating function (14) 
with y set equal to 1. The asymptotic expression (30), however, is 
rather involved and its use can be avoided as follows. Indeed, we note 
that, under the present assumptions, the instants at which overflows 
occur constitute a renewal process (i.e., the intervals between any 
pair of consecutive overflows are independent of each other and have 
the same distribution). Then using Smith's extension of a result due to 
Feller (cf. Ref. 4, pp. 296-298 and Ref. 5, pp. 30-33), we have: 

Var [S(t)] ~ fo(c) - mi 2 (cWmi 3 (c) (32) 

where ju»(c) is the nth moment of the interoverflow distribution of a 
group of c trunks. 

The expression on the right-hand side of (32) is rather easy to com- 
pute, since we have the following recurrence relations (cf. Ref. 6, p. 388) 

aMn) = 2an 2 {n) + «/x 2 (n - 1), (n = 1,2, • ■ • ) (33) 

where /n _I (c) = aEi,»(a). Hence 

Var [S(t)] ~ [fe/flW^-D+Mi'fe)] t 
Ml 8 (c) 



M2 



(c-D=2E ^--Ih^ic-l-n) 
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The second moment m (c — 1 ) can be computed either by repeated use 
of (33 ) or by means of the explicit formula 

rl (e-i)...i, 

a" 

with 

(c - 1) = 1, (c - 1)„ = (c - l)(c - 2) • • • (c - n), 

(n ^ 1). 

We note that the renewal theorem used above can be applied as long 
as the input to the system is recurrent, and the other assumptions made 
here remain the same. In these more general cases, (32) still holds, but 
the moments mi(c) and ^(c) satisfy less simple recurrence relations. 
Indeed, we have then: 

w (c) T c-i(l) -Mi(c - 1) = 

and 

j u 2 (c) 7c - 1 (l) - 2 / zi(c)[ Ml (c - 1) - 70-/(1)] - M2 (c - 1) = 

where Yn(-) is the Laplace-Stieltjes transform of the interoverfiow 
distribution of a group of n trunks and y m ' (1 ) stands for the derivative 
of7m(-)atl. 

The preceding relations follow immediately from Palm's recurrences 
(cf. Ref. 3, pp. 36-38, and Ref. 7, pp. 16-22): 

7n(s)[l - 7n-l(s) + 7n-l(s + 1)] = 7n-l(s + 1), (tl = 1,2, • ■ • ). 

The standard deviation of the call congestion computed by means 
of (29) and (31) to (33) is compared in Figs. 1-5 with simulation 
results. As may be seen from these graphs, there is good agreement 
between the theoretical and observed values. 

On each one of these charts, two additional curves are also plotted, 
namely : 

(*') {Var [S(t)]\*/ER(t) as a function of the offered load. 

This expression is an upper bound for the standard deviation of the 
call congestion. However, unless the offered load is relatively small, 
it considerably overestimates this standard deviation. 

(u) {Ei. c (a)[l - E Uc (a)]\ h /{ER(t)\ h as a function of the offered 
load. 

This quantity, referred to as the binomial approximation, is a lower 
bound for the standard deviation of the call congestion. This bound 
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underestimates the latter to such an extent, however, that it is of little 
if any value. 

In view of the agreement between the observed and theoretical vari- 
ances of the call congestion, the latter are graphed in Figs. 6-8 for 
c = 1(1)10(2) 20(5)50 and 0.1 ^ a/c ^ 10. These values pertain to the 
case t = 20. The asymptotic variance of the call congestion for any 
(sufficiently large) value of t may be obtained by multiplying the vari- 
ances of Figs. 6-8 by 20//. 

Simulation results have shown that (29) — with Var [S(t)] and 
Cov [R(t),S(t)] replaced by their respective asymptotic expressions — 
give sufficiently accurate values of the variance of the call congestion 
whenever the length of the observation period, t, is such that the ex- 
pected number of offered calls, ER (t), is about 40 or more. When ER (t) 
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BINOMIAL APPROXIMATION 

o OBSERVED STANDARD DEVIATION OF CALL CONGESTION 

(INTEGERS STAND FOR NUMBER OF HOURS IN RUN) 




6 7 8 

LOAD IN ERLANGS 



Fig. 1 — Standard deviations of call and time congestions, c = 6, t = 20. 
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Fig. 2 — Standard deviations of call and time congestions, c = 10, t — 20. 
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Fig. 3 — Standard deviations of call and time congestions, c = 20, t = 20. 
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Fig. 4 — Standard deviations nf call and time congestions, c = 30, I = 20. 
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Fig. 5 — Standard deviations of call and time congestions, c = 40, I = 20. 
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Fig. 6 — Variance of call congestion. Observation period = 20 holding times 
c = 1(1)10. 

falls below 40, (29) provides us with an upper bound which becomes 
increasingly coarse as the expected number of offered calls decreases. 



IV. RELATIVE ACCURACY OF LOSS ESTIMATES 

Various measurements can be used to estimate the probability of loss. 
The principal ones are : 

(i) The number of offered calls and the number of lost (or overflow) 
calls. The ratio of the latter to the former (i.e., the call congestion) is an 
asymptotically unbiased estimate of the probability of loss. 

(ii) The time congestion. This quantity, which is an unbiased estimate 
of the probability of loss, may be either measured exactly or estimated 
by scanning the trunks at regular intervals and observing, at each 
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scan, how many trunks are busy. The proportion of scans which indicate 
that all trunks are busy is also an unbiased estimate of the probability 
of loss. 

(Hi) The carried load (i.e., the average number of busy trunks) obtained 
either by continuous observation or by scanning at regular intervals. This 
last measurement consists in observing, at regular intervals, the number 
of busy trunks. The average of these numbers, for a given number of 
scans, is an unbiased estimate of the carried load. If L stands for the 
carried load measured either exactly or by scanning, then the demand 
rate, a, may be estimated by means of the formula 

I = d[l - E x ,c(d)\. 
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Fig. 7 — Variance of call congestion. Observation period = 20 holding times 
c = 12(2)20. 
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Fig. 8 — Variance of call congestion. Observation period = 20 holding times 
c = 25(5)50. 

E ljC (d) itself is an estimate of the probability of loss. This estimate has a 
small positive bias which tends to zero as the length of the observation 
period gets large. 

Theoretical as well as observed (simulation) values of the standard 
deviations of these estimates are plotted in Figs. 9-12. (For each load, 
the simulation results given in Fig. 12 were computed from a single 
run of 500 hours. The numerator and the denominator of each ratio 
appearing in Figs. 9 and 10 were evaluated from a single 500-hour run 
of simulated traffic.) These graphs reveal typical patterns, namely: 

(i) When the offered load, in erlangs, falls short of the number of 
trunks, the loss estimates based on continuous load measurements 
have smaller standard deviations than both the call and the time 
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Fig. 9 — Relative accuracy of grade of service estimates based on hourly 
measurements of call and time congestions — simulation results. 

congestions. In the same range, the call congestion has a larger standard 
deviation than the time congestion. 

(it) When the offered load exceeds the number of trunks, the con- 
verse situation holds; i.e., the call congestion has a smaller standard 
deviation than the time congestion, and the standard deviation of the 
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Fig. 10 — Relative accuracy of grade of service estimates based on hourly 
measurements of call and time congestions — simulation results (cont.). 
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Fig. 11 — Standard deviation of the time congestion estimated by switch- 
counting. 

latter, in turn, is exceeded by the standard deviation of loss estimates 
based on continuous carried load measurements. 

The effect of scaiming on the variances of the time congestion and of 
the loss estimates based on carried load measurements is illustrated in 
Figs. 11 and 12. 

Let us assume now that the length of the observation period is such 
that (29) closely approximates Var [S(t)/R(t)]. Under these conditions, 
the load beyond which the time congestion is less accurate (in terms 
of its variance) than the call congestion is approximately equal to the 
load a determined by the following equation: 

1 + Ei, e (a) = 2E lte (a)[c + 1 - a + atf lfC (a)]. (34) 

This condition is readily seen to be equivalent to the requirement 

Var [S(t)/R(t)] = Var Bit) (35) 

where B(t) stands for the time congestion in an observation period of 



ACCURACY OF LOSS ESTIMATES 



1161 



length t. Equations (35), (29), and (31) together with the relations 
(cf. Ref. 3 p. 131) 

ESQ) = atEB(t) 

VarS(0 - ES(t) = (at) 2 V&rB(t) 

imply (34). 

For given c, (34 ) has a unique positive root, r, which is smaller than 
c except in the case c = 1 where the root is equal to 1. Computations 
show that this root lies relatively close to c (cf. Fig. 13). 

Let B n (t) be the estimate of the time congestion obtained by switch- 
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Fig. 12 — Relative accuracy of loss estimates based on call and time conges- 
tions and on carried load measurements. 
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Fig. 13 — Root of equation (34). 

counting at the rate of n scans per observation period of length t. We shall 
now derive an explicit formula for the variance of B n (t). 

Let t be the interval separating consecutive scans, N(u) be the num- 
ber of busy trunks at time u, and 



X(u) = 



1 if N(u) = c 
if N(u) < c. 



Then 



B n (t) = n~ l EX(ir) 



and 
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EB n (t) = E uc (a) 

Varfi n (0 = n~ 2 £ E Cov [X(z»,X0V)]. 

.=i y=i 

Now let 

P( M ) = Pr[iV(«) = c |iV(0) = c]. 

The function P(- ), which is called the recovery function of the proc- 
ess N(- ), has the following expression (cf. Ref. 3, p. 85 and Ref. 8, p. 
135): 

P(u) = E ilC (a) - E — II fl — ) 

Pi Wj ttj \ Wj - WiJ 

where, as before, w< , t = 1, • • • , c, are the c roots of <r c (w + 1 ). Since 

Cov [X(u + v),X(v)] = tf 1>e (a)P(w) - ft,/ (a) 
we have (cf. Ref. 8, pp. 136-138) 



n n 



(36) 



Var B n (t) = n" 2 ft, e (a) £ £ P(l * " 3 I " A., to 

»-=i y=i 

n c |fc|u)jT 

= - n" 2 ft., (a) £ (* ~ I * I) E V- 

•n6 — M 

= n" 1 ft„(a) £ I 1 II (l — )} 

If we let n tend to infinity in this formula, we obtain, in the limit, the 
variance of the time congestion, B(t), for continuous observation (meas- 
urement): 

+ i|;»r'(i-e"og(i-^)}. 

This last formula was first obtained, in a slightly different form, by 
Kosten, Manning and Garwood (cf. Ref. 1). 
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